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Abstract 

We consider a possible mechanism of thermalization of nucleons in relativistic heavy-ion collisions. Our 
model belongs, to a certain degree, to the transport ones; we investigate the evolution of the system created 
in nucleus-nucleus collision, but we parametrize this development by the number of collisions of every particle 
during evolution rather than by the time variable. We based on the assumption that the nucleon momentum 
transfer after several nucleon-nucleon (-hadron) elastic and inelastic collisions becomes a random quantity 
driven by a proper distribution. This randomization results in a smearing of the nucleon momenta about their 
initial values and, as a consequence, in their partial isotropization and thermalization. The trial evaluation 
is made in the framework of a toy model. We show that the proposed scheme can be used for extraction of 
the physical information from experimental data on nucleon rapidity distribution. 

1 Introduction 

The problem of isotropization and thermalization in the course of collisions between heavy relativistic ions 
attracts much attention, because the application of thermodynamic models is one of the basic phenomenological 
approaches to the description of experimental data. Moreover, the assumption about a local thermodynamic 
equilibrium, along with other factors, is successfully used in various domains of high-energy physics. Meanwhile, 
many of the questions concerning this issue remain open for discussion [l]-[8]. 

The main goal of the investigations of the collisions of relativistic nuclei is the extraction of a physical 
information about nuclear matter and its constituents. It is a matter of fact that we can get know more about 
quarks and gluons (constituents) just under extreme conditions, i.e. at high densities and temperatures. During 
last two decades, one of the celebrated tools on the way of the theoretical investigations of this subject was 
relativistic hydrodynamics which started to be applied to elementary particle physics by the famous Landau's 
paper [9j. 

Applying relativistic hydrodynamics, one can partially describe experimental data and get know that the 
matter created in relativistic nucleus-nucleus collisions can be regarded on some stage of evolution as a continuous 
one, i.e. as a liquid. Moreover, as was discovered in Brookheaven National Laboratory, it can be regarded 
even as a perfect fluid [TUIIII] which consistent with a description of the created quark-gluon plasma (QGP). 
Unfortunately, the important physical information is hidden in sophisticated numerical codes which solve the 
Euler hydrodynamic equations of motion. 

In the present paper, wc propose an approach to description of relativistic heavy-ion collisions which allows 
one to extract the physical information from experimental data on the basis of a transparent analytical model. 



2 The Model 

Our model is aimed at the extraction of the physical information from the nucleon spectra for such collision 
energies when the number of created nucleon-antinucleon pairs is much less than the number of net nucleons. 
This means that the model can be applied at AGS and low SPS energies. 

1. We separate all nucleons in the final state, i.e. after freeze-out, into two groups in accordance with their 
origination: a) the first group consists of net nucleons that went through just hadron reactions; b) the second 
group includes nucleons which were created in the collective processes, for instance, during hadronization of the 
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Figure 1: Transformation of the initial nucleon 
momentum, fco, as a result of M collisions; Q = 
X]m=i 9'™'' ^^'^ total momentum transfer after 
M effective collisions, where q'-™^ is the momentum 
transfer in the m-th collision. 



Figure 2: Example of inelastic collision of two nu- 
cleons with creation of two 7r-mesons, TV + — > 
N + N + n + TT. 



QGP. In accordance with this, the nucleon momentum spectrum can be represented as a sum of two different 
contributions: 



dN fdN 



- > , ■ > -1 • (1) 

d^p \d^pj^^^^^^ \d^P J QGP 

In turn, the total number of registered nucleons equals A^totai = -/Vhadron + ^qgp- If this separation can be 
done, then we can define the "nucleon power" of the created QGP as -Pqqp = -/Vqcp/Motai- In our further 
investigation, we will deal mainly with the nucleons from the first group which come to the final state after a 
chain of sequential nucleon-nuclcon (-hadron) elastic and inelastic collisions. 

2. The collision number for every nucleon (hadron) is finite because the lifetime of the fireball is limited. To 
determine the maximal number of collisions, M,„ax, in a particular experiment we use the results of UrQMD 
simulations [121 [T3] . 

3. Because the colliding nuclei are the spatially restricted many-nuclcon systems, the different nucleons 
experience different collision numbers: it is intuitively clear that the collision histories of the inner and surface 
nucleons will be different. That is why, we partition all amount of nucleons of the first group into different 
ensembles in accordance with a number of collisions before freeze-out. Then the nucleons from every ensemble 
give their own contribution to the total nucleon spectrum. If we denote the number of particles in a particular 
ensemble where the particles experienced M collisions by C(M), then, in correspondence to Eq. ([1]), we can 
write the total nucleon spectrum as 

= V C{M)Dm{p) + CthcrmZ?thcrm(p) , (2) 

where Dm{p) is the spectrum (normalized to unity) of the particles in the M-th ensemble. The last term on the 
r.h.s. of ([2]) corresponds to the possible contribution from the totally thermalized source which we associate with 
QGP. Here, Ctherm is the number of nucleons from the second group which are created during the hadronization, 
and £'tiicrm(p) is the thermal distribution normalized to unity. 

Consider successive variations of the momentum of the n-th nucleon from nucleus A which moves along the 
collision axis from left to right toward nucleus B. Every jn-tli collision induces the momentum transfer, gn™\ 
for the n-th nucleon. So that, after M collisions, the nucleon acquires the momentum A;„: 

fco ^ feo + gi^^ ^ k^ + q^n^+q^n^ ^ ■•• ^ fco + Q„ = fc™, (3) 

where (5„ = ^ <7„ is the total momentum transfer finally obtained by the n-tli nucleon (see Fig.[T]). 

rn— 1 

Let US assume for the moment that the elastic scattering gives the main contribution to the two-nucleon 
collision amplitude. The initial momentum of every nucleon in nucleus A \& ka ^ = (0,0, /cqz); while the 
initial momentum of every nucleon in nucleus B is fct, = = (0,0, —k^z)- The energy and momentum are 
conserved in every separate collision of two particles, uj{ka) +uj{kt) = to{Pa) + '^(Pb), ka + ki, = p^+p^^ where 
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and are the momenta of the particles after the colhsion. We assume that the particles are on the mass 
shell, so that uj(k) = \J + (the system of units ^ = c = 1 is adopted). Determining the six unknown 
quantities, and p;,, from four equations is straightforward, but two quantities, e.g. {Pa)x and {Ph)x, remain 
uncertain and can be considered as such which accept random values driven by the scattering probability. So, 
after the first collision, one component of the particle momentum becomes random. After the third collision of 
the particle, each component of the momentum becomes random, hence the particle momentum, for instance 
p^, becomes completely random. As a result, we can regard the momentum transfer, q = p^ ka, as a random 
quantity as well. If we follow the nucleon sequential elastic scattering from the first collision to the last one, we 
would see the full randomization of the particle momentum after each three sequential acts of scattering. So, 
we can look at the process of randomization from the equivalent point of view: starting from the definite initial 
momentum fco, the n-th particle gains completely random momentum transfer qiP, Qr? ■, ■ • • , Qri"^^ after each 
three sequential collisions in the chain of 3M collisions, see ([3]). 

At the same time, the nucleon momentum transfer undergoes the even faster randomization in the inelastic 
collisions. Indeed, let us consider, for example, the process of creation of N^^ pions in the nuclcon-nucleon 
reaction: N + N ^ N + N + Nj^ ■ tt. The process leads to randomization of [3(2 + iV^) — 4] degrees of freedom, or 
dp degrees of freedom per particle becomes random, where dp = 3 — 4/(2 + A^^) and we assume that all particles 
are on the mass shell. For instance, if Nt^ — 2 (sec the physical diagram of the collision depicted in Fig. [2|) , then 
two components of the momentum of every particle after a reaction becomes random, i.e. dp = 2 in this process. 
If ^ 1, the number of the random degrees of freedom per particle achieves its maximal value dp = 3. Thus, 
we see that, in the inelastic collisions, the randomization of the nucleon momentum transfer attributed to one 
physical diagram goes faster than that in the elastic scattering. 

To estimate the effective number of collisions, M, we analyze all nucleon collisions (physical diagrams), 
Ncoii, and obtain a total sum of the random degrees of freedom, dtotj gained by the particle during all reactions 
before freeze-out, i.e. we need to know dtot = J2f=i^ ■ Then, the effective number of collisions is determined 
as M = fitot/3. So, we assume a randomization of the momentum transfer, q, after every effective collision; the 
number of effective collisions, M, is determined by the conditions of a particular experiment. 



2.1 Many-particle distribution function for M-th ensemble 

We would like to determine the density distribution function in the momentum space, f2N{Etot] ki, . . . , fe2jv)i 
which describes 2A" nucleons in the final state, after M effective collisions which are experienced by every 
nucleon before freeze-out (we name this group of nucleons as the Af-th ensemble). All consideration is carried 
out in the c.m.s. of two identical colliding nuclei, hence the initial total momentum of the system A + B is equal 
to zero. Let us write down the density distribution function in the form 



/2Jv(£^tot; fcl, • ■ • , fc2Ar) 
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f2N{Etot\ki, . . . , fc27v) 



(4) 



where we normalize the density distribution function in such a way that simultaneously determines the density 
of states in the system: 



^2N{Etot) 



dki . . . dk2N /2Jv(i?tot; fci, . . . , fc2Jv) . 



(5) 



The measure of integration in the single-particle phase space looks like (in units of h) dkn — ^■^^ • 

The unnormalized distribution function /2w(£'tot; fci, ■ • ■ , fc2Ar) is determined in a two-fold way: first, we 
follow all collisions of a particular nucleon by integration with respect to all nucleon random momentum transfer, 
and, second, we fix the total energy of the 2A-nucleon system after freeze-out in a microcanonical-like way: 
Etot = J2n=i^i^ri)- Then, it reads 



2N 



f2N{Etot;ki, . . . ,fe2Ar) = si Etot - ^ w(fcn) 



dPi dP2N TT 

n=l 



V 



V 



M 



i2nfS^kr.~ko~Y. 1 



(m) 



2N 

n 



M 



(6) 



i2nf (5^ [K.,+ko~Y. 1 
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where V is the volume of the system in the coordinate space. Here, we assume the independence of the sequential 
scatterings, which results in the independence of the particular momentum transfer (see for details [14j). Hence, 
the element of the probability to accept a particular chain of the momentum transfer by the n-th particle in 
the scries of M collisions reads 



M 



m— 1 



(m) 



(m) 



(27r)3 



(7) 



where the distribution of the momentum transfer in the m-th collision is characterized by the presence of the 
form- factor J,„(<7) with 



Jvi{q) 



1 . 



(8) 



In what follows, for the sake of simplicity, we assume the independence of Jm{q) on the collision number, 
i.e. Jra (g) — * {Jm (?) )coiiisions ~ (*?) ■ Heucc, wc adopt an approximation where just one form-factor J{q) 
characterizes a distribution of the momentum-transfer in a series of collisions which are experienced by a nuclcon 
during its traveling through the fireball. 

With making use of this approximation, we make a step in the evaluation of the 2A^-distribution function 
. If one represents (5- functions in ([5]) in terms of the Fourier integrals with respect to the auxiliary variables 
a„, hn and then with making allowance for definition of the integration measure ([7]), the unnormalized density 
distribution ^ can be written down in the form 



2N 



N 



f2N{Etot] ki 



2N) 



n=l 



27V 

n 



ri=JV+l 



V 



-ian ■ {hn—ko ) 



M 



-ibn ■ {kn-\-ko) 



M 



(9) 



where J{r) = J J{q) e*'^ '' is the Fourier component of the momentum-transfer form-factor which corre- 
sponds to the single nucleon scattering. 

We introduce the "multi-scattering form-factor" (M-scattcring form-factor) 



hiiQ) 



1 

V 



d'^re' 



[J{r)] 



M 



Then we can rewrite the unnormalized distribution function 

/ 2N \ N 



2N 



f2N {Etot ; fei 



2N) 



fco) 



(10) 



(11) 



n=l 



n=l 



n'=N+l 



To obtain the partition function, we make the Laplace transformation of the density-of-states function ([5]) 
with respect to the variable Etot (the total energy of the nuclcon system after freeze-out). As a result, we obtain 
the partition function of the canonical ensemble 



oo 

Z2nW= J dEtote-P''^°^n2N{Etot). 



(12) 



If wc define the unnormalized 2A^-particle distribution function of the canonical ensemble as 

oo 

F2Ar(/9;fei,...,fc27v) j dEtote~'^^'°'f2N{Etouki,...,k2N), 



we can write the partition function as the integral 
Z2n{P) = / dki . . .dk2N^2N{f3;ki, . . . ,k2N) ■ 



(13) 



(14) 
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It is obvious that the partition function Z2N iP) plays the role of the normalization constant with respect to the 
2Af-particle distribution function F2Ar defined in p3p . Next, we determine the 2A^-particle distribution function 
in the canonical ensemble as 

F2Ar (/?; fcl, . . . , k2N) = „ ^2n{I3; fcl, . . . , k2N) ■ (15) 

Taking Eq. (fTTj) into account, for the distribution functions F2Ar(/3; fei, . . . ,k2N) from (fTS]) . we obtain 

N 2N 

¥2n{(3; fcl, ... , k2N) = n [e-^"^'=")/M(fe„ - feo)] n [e-^^^""'^ hi{k.n' + feo)] . (16) 

n=l n'=N+l 

So, we can write the distribution function which characterizes the 2iV-nucleon system where each particle 
experiences M effective collisions (A/-th ensemble) in a factorized form 

N 2N 

¥2N{P,ki,...,k2N)^Y[fa{kn) J] Z''^'^")' (1^) 

n=l n=N+l 

where 

fa{b){k) = —- e-^-^'^) lM{k T fco) with z,(,)(/3) = V / —3 e-'^-C") /^(fc T feo) , (18) 

are the single-particle distribution function and the single-particle partition function, respectively, attributed 
to nucleus "A" for subindex a or to nucleus "B" for subindex b. 

2.1.1 Saddle-point approximation 

For large enough collision numbers M, we can calculate the integral in ()10p within the saddle-point method. 
To do this, we first represent the integral in the form Im{Q) ~ y J d^T exp [~iQ ■ r + M\nJ{r)] . Then we 
expand the exponent In J(r) into a series at the point r = up to the second order. We take into account that 
J(r)|^_P = 1 (see Eq. ([8|)). If the form-factor depends on the modulus of the momentum transfer J(|g|), the 
expansion of the logarithm In J(r) looks like 



lnJ(r)« -i(g2)r2, where {q^) ^ J jQq\) (19) 



with q = \q\ and r = \r\. Here, we take into account that 5,;J(r)|^ ^ = 0, which is valid when J{—q) = J{q)- 

Finally, we obtain integral (jlOp in the approximate form 



where Q = \Q\. 

So, the M-scattering form- factor Im{Q) for a large number of the effective collisions M can be approximately 
represented as the normal distribution oc exp (— (5^/2(t^) with the variance a ~ ^ M{q^) /i , which is obviously 
a reflection of the central limit theorem. 

Following approximation ([20]) . the single-particle distribution functions (fT8|) take the form 

fa{b){k) = — e ^^n^ with ZM /3 = / e ^I7(^ , 21 

where we skip the common factors in fa(b){k) and in ZM^fi), respectively. In the limit case M — > 00, the 
dependence on the initial momentum ±fco is washed out, and both single-particle distributions fa{k) and fb{k) 
take the same "thermal" limit: 

faib)ik) - /thc„„(fc) e-^-('=) , zthe™(/3) = / e-^-e^) . (22) 
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In the non-relativistic case, i.e. when Lo{k) = /2m + m, one can rewrite distributions (|2ip as the Jiittner 
ones [15] (see [T6] ) 



fa(b){k) 



exp 



(fc =F mUh)^ 



with zm{I3) ~ 



Here, we define the collective (hydrodynamical) velocity as 



(fik 
(2^ 



exp 



(k — mt^h)^ 



(23) 



k 



t3 

/3coll(M) , 





m 



with 



3m 

M{q^) 



and the effective temperature as 



Tcs — 



1 



/3 + /3coii(M) 



(24) 



(25) 



The quantity PcoW can be put in correspondence to the "collision" temperature Tcoii = l//3coii = f-^-^(^(9))i 
where (w(q)) = {q'^)/2m is the mean energy transfer in one nuclcon collision. We see that every collision 
ensemble has its own effective temperature and own collective velocity. Indeed, with increase in the number 
of collisions M, the effective temperature increases and the collective velocity is going down. These quantities 
have the following limits: 



(M = 0) 



kg 



> 



> 



< Teff < T 



(M oo) 



(26) 



where w(fco) = ym^ + fcoi ~ V/^j ^^^i we marked conventionally the nucleons which do not take part in 
any collisions (spectators) by M = 0. Here, the left limit corresponds to the case M ~ 0, and the right limit 
corresponds to AI —> oo. Hence, the every step during the increase in the number of collisions results in the 
redistibution of the energy accumulated initially in the longitudinal movement, i.e. Vh = fco/w(fco) and Teg = 
before first collision. The energy passes partially to the transverse degrees of freedom, which increases the 
isotropization and the temperature of the system and decreases, in turn, the collective velocity. When the 
number of collisions is large enough, i.e. when M oo, we come to the limits: ^ and Tcff ^ T. These 
conclusions are valid, of course, for relativistic energies as well. 

Meanwhile, when the number of collisions is finite we see that, just for the non-relativistic dispersion law, 
the obtained distribution (pij) can be regarded as a locally thermal one (the Jiittner distribution function). 



2.1.2 Single-particle spectrum for M-th ensemble 

We can construct a "two-source" single-particle distribution function f{ka,kt,) — fa{ka)fb{kb) 
averaging of the two-source random quantity W{ka, kb) gives 



(W) 



d^ka d^kb 



(27r)3 (27r)3 



Wika,kf,)fika,ki,). 



Then the 



(27) 



If we take W{ka, kb) = ^ [S^{p — ka) + 5^{p — kb)\ , we obtain, after the averaging, the single-particle spectrum 
which is attributed to the system where all particles experienced M collisions before freeze-out: 



M 



h-iip-ko) hi^p + ka) 



(28) 



where Im{Q) is defined in ([TO]). Hence, with accounting for notations (flS)) . the distribution functions Dm{p) 
can be represented also as 



DM{p)^\[fa{p)+h{p)] ■ 

For large collision numbers M, the distribution function D]\j{p) reads 



Dm{p) 



2zm{(3) 



3(p-fco)^ 3(p + fc0)- 



with ZMiP) 



d^k - 

re 



(27r)3 



(29) 



(30) 
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It is evidently seen that the spectrum has two items which can be attributed to the first and to the second 
colliding nuclei, respectively. In accordance with this structure, it can be named a "two-source single-particle 
spectrum" . It is clear also that this picture is just the effective one: after a nucleon-nucleon collision, wc cannot 
distinguish nucleons, and thus wc cannot say anything about their origination from a particular nucleus. 

In what follows, wc use the distributions Dm{p) from (j28|) to describe a nuclcon spectrum arising in rcla- 
tivistic nucleus- nucleus collisions. 



2.2 Nucleon Rapidity Distribution and Transverse Spectrum 

To obtain the transverse mass and rapidity distributions, we pass to new variables {px,Py,Pz) — ^ 

m± = {rn^ + p^Y^'^, p\ ~ Py^ tanhy = ^ , then (Pp = dcjyujp m± dm± dy, where (f> is the azimuth 

angle. In accordance with the double differential cross-section reads (for the sake of simplicity, we consider 
here the most central collisions which possess the azimuth symmetry) 



= ZTT m± cosh y 



m±dm±dy 



Mm 



C{M) DM{m±,y) + Cthcrm Allcrm("^J., y) 



M=l 

We can define the distribution functions in the rapidity space as 



27r cosh y 

2:therm(/3) (3^ COSli' y' 



with 



X = mPcoshy and 2:thorm(/3) = 47ri^ K2{ml3). 



(31) 



$M (y) = 27r cosh y / dmj_m\ DM{mj_,y) with / dy$A/(y) = l. (32) 

J m J 

Then the rapidity distribution looks like 

— ^ ^ C(Af) (2/)+ Cthcrm Cthcrm (y), (33) 

^ M=l 

where 



*thcrm(2/) = 7 7^ ,3 ix^ + 2x + 2) (34) 



3 Toy model 

We would like to obtain the explicit results in the framework of the proposed approach. We note that the crucial 
quantity of our approach is the form- factor J{q) which can be considered as the probability density (see ([5])) to 
find the nuclcon momentum transfer q in a particular nucleon collision. Then the mean value of the quantity 
which depends on the momentum transfer, for instance /(q), is obtained as (/) = / d^q/ {211)^ f{q) J{q). To 
estimate the behavior of the form- factor J(q), we made evaluation of the mean values of g^, and q1 in the 
framework of the UrQMD [I2l[l3j. The results of the evaluation are depicted in Fig. [3l It is evident that, for 
the central collisions, the equality (g^) = (g^) should be valid; we see this, indeed, in Fig. [3l {q\) « 2((7^). It 
is seen that, for AGS energies, the mean values (g^) a-nd (g^) approximately equal one another or of the same 
order. We take this as a basis to formulate a model: let the form-factor J{q) be chosen as a homogeneous 
distribution of the momentum transfer in the sphere of finite radius gmax (see Fig. |4]), 

(2tt)^ 4 r d^q 

J{q) = d{qme.K -\q\), where Vg ^ -nq^^^ and J J{q) = 1 . (35) 

As a matter of fact, an application of this toy model can be regarded just as an approximate description of the 
nucleon distributions at AGS energies. 

We calculated the Fourier transform of form-factor ([35|l in the explicit way and obtained 

3 

J{x) — — ( sin X — X cos x ) , (36) 
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Figure 3: The mean value of the nucleon momentum 
transfer squared as a function of the coUision energy 
obtained in the framework of the UrQMD [121 [13] sim- 
ulations for the most central collisions at AGS ener- 
gies. 



Figure 4: Nucleon momentum transformation as a re- 
sult of the two-particle collision; k is the initial mo- 
mentum, q is the momentum transfer, fc + q is the 
momentum after the reaction, (/max is the maximally 
allowed modulus of the momentum transfer. 



where x = rgmax and r = \r\. So, the multiscattcring form-factor Im{Q); see (ITUl) . reads now as 

"3 1 

— (sinx — X co&x) 

X"^ 



hi{Q) = 777: / drr s\n(Q'i 



(37) 



where Q = \Q\, x = rq^^ax- 

With the help of these functions, we can evaluate the partial distribution Dm{p) from ([28]) as 



1 



2z(/3) 



-/3m_i_ cosh y 



[lM{\p-ko\) + lMi\p + ko\)] 



where = m± coshy, = m± sinh y, |p =p feol = v^m^ — + {m± sinhy ^ fcoz)^ and 



z(/3) = V 



(27r)3 



Iai{\p - ko\ 



(38) 



(39) 



Next, we use these partial distribution functions directly in ([3T|) to evaluate the double differential (transverse 
mass) spectrum and the functions ^i\j{y) from (|32p which will be used then as the partial rapidity distributions 
in ([351). 

In the framework of the toy model, we obtain that the parameter of the model, qmax, can be related to the 
mean value of the momentum transfer squared 



/ 2\ ^ ' 



2 

max • 



(40) 



Taking into account {q'^) evaluated in UrQMD [121 \IM AGS energies (see Fig. [3]), we obtain the following 
estimation from (gOl): qmax ~ 800 McV. 

In the framework of the toy model for large collision numbers M, the distribution functions DM(m±,y) 
reads 



-I3m± cosh y 



2z{P) 



5(p + fco) 



with z(/?) 



^e-^-(^)-l^ (41) 
(2^)3" ' ^^'> 



where |p — feo| and |p -I- feo| are defined in the same way as in 
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3.1 Extraction of the physical information from experimental data 



By making use of the rapidity distribution (j33p , we fit the experimental data on the rapidity distribution of net 
protons which were measured at the AGS (E802 Collaboration) [T7]- The proton data is remarkable in that 
sense that we know exactly the initial momentum, fcoz, of every nucleon. 

We note that, in the case of a small number of experimental points iVcxp, the set of functions $m(2/) is 
overcomplete (A'oxp < Afmax + !)• To choose a unique configuration of the variable parameters C(M), we use 
the maximum entropy method |181 119| . Namely, we define the information entropy 



cr = - ^ C{M) lnC(A/), 



(42) 



M=l 



where C{M) = C{M)/Np, Np is the total number of protons, and, for the unification of notations, we adopt 
C(Mniax + 1) = C'thcrm- Our goal is to find the maximum of the information entropy a. Meanwhile, entropy 
(j42l) is supplemented by (A^oxp + 1) constraints: 



M„ 



(cxp) _ 



K+1 



E C{M)<PM{y^).. 



A/=l 



where /I'^'^p' is an experimental data value of the distribution dN/dy taken at the rapidity point yi, i 



1,2, 



, TVexp and 



^+i{y) = *&therm(y)- Usiug thc mcthod of the Lagrangian multipliers, we reformulate 



the problem in the following way: it is necessary to find the maximum of the expression 



L 



A/max + l 

E 

A/=l 



C{M) h-iC{M)+fi 



M„ 



E cw-i 



A/=l 



Wcxp 

E- 

i=l 



r(cxp) 



A/„ax + l 

E 

A/=l 



C{M) $A/(yO 



where fi and A,; are the {Ncxp + 1) Lagrangian multipliers. Derivatives of this expression with respect to 
all unknown coefficients give a set of (Afmax + 1 + -^oxp + 1) equations which we solve using the variation 
method: first, we exclude and all C(M) from the system of equations and left with a reduced system of 



iVexp transcendental equations for the unknown Lagrangian multipliers Ai, where i = 1,2,. 



The latter 



reduced system of equations is solved by variations of the Ai. So, we seek out the minimum of the expression 



Wo: 



x'(Ai,A2 



>) = E 



r(cxp) 



-N, 



A/max + 1 

. E 

A/=l 



C{M) <^>M{y^) 



(43) 



where C{M) = Xm/ (Ea/=T'''^ ^a/) with Xm = cxp Np J^fj^l^ Xi'^M{y^) ■ Actually, Eq. (O is nothing more 

as a total x^i with a help of which we find a theoretical curve of the closest fit to dN/dy experimental data. 
At the same time, it is easy to rewrite (j43|) with accounting for normalization with respect to the experimental 
error bars at every rapidity point yi. 

The fit was carried out with a help of the program MINUIT, and the variable parameters are the Lagrangian 
multipliers A^. The coefficients C(M) and Ctherm are evaluated then through the obtained multipliers A^. 
Thc slope parameter f3 was first extracted from the double differential yield for protons with the use of the 
thermal distribution. All evaluations are carried out in the c.m.s. of colliding nuclei. The obtained theoretical 
curves are depicted in Fig. [S] The solid thick curve represents the result of the fit, total = 13-3. The 
broken curves marked by the numbers M represent the partial contributions from every ensemble, C{M) (pM{y)- 
The thermal contribution is represented by the central bell-like dashed curve. The obtained parameters for 
T = 1//3 = 280 MeV are shown in Table 1. In the form of histograms, these coefficients are depicted in Fig. [5] 

Table 1 





C(2) 


C(3) 


C(4) 


C(5) 


C(6) 


C(7) 


C(8) 


C(9) 


C(10) 


C(ll) 


C(12) 


C(13) 


C'thcrm 


1.3 


23.2 


4.7 


5.3 


6.6 


8.2 


9.5 


10.6 


11.6 


12.3 


12.9 


13.5 


13.9 


16.6 



So, each coefficient C(M) tells us how many protons undergo M effective collisions or which is the population of 
each collision ensemble. For instance, the ensemble of protons which participated just in two effective collisions, 
M = 2, consists of 23 protons, i.e. C(2) = 23.2. It is worth to note that this ensemble is maximally populated 
in comparison to others. What is very important, we learn from this expansion that Cthorm = 16.6. This 
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Figure 5: The result of the fit (thick soUd curve) 
to experimental data [17j on rapidity distribution. 
Broken curves marked by numbers M represent the 
partial contributions from every collision ensemble, 
C{M) ^M{y), the thermal contribution is represented 
by the central Gaussian-like dashed curve . 



AGS (£802): Au+Au, p,^^=11.6A GeV/c 




I < I < I < I < I < I < I " I i-r' I ' I ' I < I < I < I < I < I I 
012345678 9 10 11 1213141516 



M 



Figure 6: The values of the coefficients C{M) which 
were obtained as the result of a fit of the dN/dy data 
[IT] . The histogram of the coefficient Cthcrm is putted 
on the place M = 15 just for a convenience of com- 
parison. 



means that approximately seventeen protons come from a thermal source, what makes up 11% (eleven percent) 
of all participated protons. The maximal number of collisions, M„iax = 13 was determined from the UrQMD 
evaluation [121 113| as the maximum number of effective proton collisions which is obtained by averaging over 
the number of Monte-Carlo events. Actually, the histograms depicted in Fig. [H] represent a tomography picture 
with respect to the number of collisions. 

The partial distribution functions in the rapidity space, $m(2/) and $therm(y) (see ((32)) and ([34)) ). are depicted 
in Fig. [T] One can say that the total rapidity distribution dN/ dy from (jSSp is nothing more as an expansion 
over the set of these A/max + 1 functions. In Fig. [71 we depicted also, as the scatter curves, the functions <&A/(y) 
evaluated in the framework of the saddle-point approximation (see ([2T|)). Starting from M ~ 4, these curves, 
as seen, differ very slightly from the distribution evaluated in a rigorous way. 

So, we sec that the top of partial rapidity distributions which corresponds to the rapidity point j/top is shifted 
to the central rapidity region with increase in the number of effective collisions M (in Fig. [71 the number AI marks 
every particular distribution) . This means that the collective velocity ~ tanh ?/top of the distribution <I>j\/ (y) 
decreases with increase in M . Evidently, the energy which was initially accumulated in the longitudinal degrees 
of freedom is converted to the particle creation and to the excitation of the transverse degrees of freedom with 
decrease in the longitudinal velocity Vh, i.e. after every collision. The latter results in the collision- by-collision 
increase of the isotropization and the effective temperature. These two conclusions are in full correspondence 
with the analysis made in Section r2. 1.1 1 

By making use of the obtained coefhcients Cm and Ctherm, which are listed in Table 1 (see histograms of the 
coefficients in Fig. [6]), we evaluated the double differential cross-section in accordance with ([3T|). The results of 
this evaluation for different rapidity windows is depicted in Fig. [51 A very good description of the experimental 
data [l7] for all rapidity windows is seen. Hence, our model satisfactory describes the m^-spectra. 

We assume that the thermal source has absolutely different nature of origination, i.e. it cannot be created 
just due to the hadron reactions of nucleons which result in the randomization and the subsequent isotropization 
of the nucleon momentum. The thermal source can emerge as a result of the appearance of many other (not 
hadron) degrees of freedom. We know just one candidate to this role, it is the quark-gluon plasma, for instance, 
its creation can occur in collisions of nucleons at high energies in the presence of a dense medium. A many- 
parton system, which emerges in the collision, consists of a large number of gluons and quarks. All momenta 
of quarks and gluons can be regarded from the very beginning as random ones, and the thermalization of the 
system occurs during a time span Tthorm = 0.6 fm/c [20j . Hence, the protons which come from the thermal 
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source indicate the presence of the QGP in the fireball, and we can determine a nucleon power of the QGP, -Pqgp, 
by the ratio of the number of protons coming from the thermal source, Cthorm, and total number of participated 
protons: Pqgp = Cthovm/Np. For instance, it turns out that the nucleon power of the QGP equals eleven percent, 
Pqgp « 11%, in Au + Au collisions |T7] (0-3% centrality). 
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Figure 7: The curves marked by numbers M rep- 
resent the functions <I>a/(2/) (see ([32]) ). the thermal 
distribution $thcrm(y) is represented by the central 
Gaussian-likc curve. The scatter curves represent the 
same $m(2/) evaluated in the saddle-point approxima- 
tion. 



Figure 8: Solid curves represent the evaluation of the 
TO^-spectra obtained in accordance with Eq. (|3ip . 
where we use the same values of the coefficients C{M) 
which were obtained as the result of a fit of the dN / dy 
data (see Fig. [5]). Experimental points arc from [17] . 



4 Summary and discussion 

The starting assumption of our approach is given in Eq. ([T]) : we propose to separate the registered nucleons by 
their origination into two groups. The first group consists of the net nucleons after the multiple rescattcring 
and the second group consists of the nucleons created by the QGP. We consider the latter as a source of the 
fully thcrmalizcd nucleons which is in the rest in the c.m.s. of colliding nuclei. Our model is applied at AGS 
and low SPS energies, where the creation of nucleon-antinucleon pairs is comparatively small. 

The second assumption consists in full independence of the momentum transfer in the chain of sequential 
collisions of a nucleon: 

M 
m— 1 

where J (q'-^-' , . . . , q'-*^'') is the probability density to find the set q*^^^ , • ■ • , q*-^^-* of the momentum transfer during 
M collisions of the nucleon in the fireball; Jm{q) is the probability density to find the momentum transfer q 
in the m-th collision. Additionally, we assume that the scattering conditions do not depend on the collision 
number: 

Jm{q) ^ (-/m (g) )comsio„s = (?) ' 

The essence of the third assumption is a partition of all nucleons, which took part in rescattcring processes, 
into collisions ensembles in accordance with the effective collision number M, see Eq. ([T|). The physical basis of 
this assumption is a spatial finiteness of the many-nucleon system. 

So, we parametrize the time axis by the number of particle's effective collisions M = (i')t, where (v) is 
the mean frequency of collisions and t is the time interval. If we considers the dependence of the distribution 
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function ([^5]) or ([50]) on the variable M, we can regard this distribution as a nonequihbrium one {M oc time). 
As we see from ((30)) . the increase of M effectively mimics an approach to the Boltzmann distribution (for the 
sake of simplicity, we consider no special statistics). If the number of collisions M is fixed, this means that the 
process of thermalization stopped at the time moment when the particles had experienced just M collisions 
and, at this moment, they were frozen out; all this results in a partial thermalization of the subsystem which 
we name as the M-th collision ensemble. So, the level of thermalization and isotropization depends on the 
number of effective collisions, M, which is determined, first, by the lifetime of the system or, more precisely, by 
the number of physical collisions of every particle and, second, by the level of randomization of the momentum 
transfer in every physical collision. 

Actually, in the framework of the proposed model, we obtain a physical picture which looks like a discrete 
"fireball model", where every fireball can be associated with a particular collision ensemble. In accordance 
with the fireball model [21| . the rapidity axis is populated with thermalized fireballs following a distribution 
PiUh) which is taken as the Gaussian one. (Here, the rapidity yh determines the collective (hydrodynamical) 
velocity Vh = tanhj/h of the particles attributed to the particular fireball.) The essential difference which occurs 
in our model, except the discretization of the fireball set, is that "our" fireballs (collision ensembles) are not 
yet the fully thermalized systems. The degree of thermalization of every collision ensemble is frozen on some 
stage of evolution toward a full thermalization (the collective velocity of the particles which belong to the 
particular ensemble can be evaluated approximately as v\i « tanhj/top, where ytop determines the top of the 
^Aiiy) distribution in the rapidity space, see Fig.[7|). 

Properly, the results of the proposed model crucially depend on the single-particle form-factor J{q) which is 
nothing more as the probability density to find a particular value q of the momentum transfer in one collision. 
As a trial evaluation (toy model), we take the form- factor as a homogeneous restricted distribution of the 

momentum transfer, J{q) = 0{qmax - \q\), where Vq = ^gfnax (see Fig. |4]). 

The maximum number of collisions (reactions) Mmax is assumed to be finite and determined by the nuclear 
number A, initial energy, and centrality. With the help of the UrQMD transport model [THUS], it was found that, 
under AGS (Au-|-Au, 11. 6A GeV/c, 0-3% centrality) conditions [17], Mmax = 13. Using the thermal distribution, 
we extract the slope parameter from the experimental data on the proton mj^-spectra, T = 280 MeV. 

We made fit of the experimental data [17] on the rapidity distribution of net protons and obtained the 
collection of coefficients C{M) (see Table 1 and Fig. ^ which are nothing more as the absolute number of 
protons in every collision ensemble, i.e. Np = X]a/=T C!{M) + Cthorm, where Np is the total number of protons. 

The knowledge of the number of protons, Cthorm, which come from the QGP gives us a possibility to evaluate 
the "nuclcon power" of the QGP, Pqqp, created in a particular experiment on the nucleus- nucleus collision. 

We find that, under the AGS conditions [T7| (a centrality of 0-3%), P^^p « 11%. So, in the framework of the 
proposed criterion, it could be claimed that the QGP (as a nucleon source) was created not only at SPS energies 
|22j . but it was also created in the central collisions at AGS energies. 

Meanwhile, if we do not use the thermal contribution in the partial expansion (|33p for the description of 
the particular experimental data [T7], we obtain as well a good description (the same x^) of both the rapidity 
distribution and TO_L-spectra. So, we cannot resolve unambiguously the problem about the presence of the 
thermal source. In fact, to overcome the problem, we need a more detailed experimental information for the 
central rapidity region. 

Solid curves in Fig. [5] represent the results of calculations of the mj_-spectra obtained in accordance with 
Eq. (PT|) . where we use the same values of the coefficients C{M) (see Table 1) which were obtained as the result 
of a fit of the dN/dy data. It is seen from the results depicted in Figs. [5] and [7] that the main contribution to 
the distribution for small rapidity values (in the laboratory system) is given by collision ensembles with small 
M. For instance, for the particular experiment [IT], the partial distribution for M = 2 in Eq. ([33)1 determines 
the distribution in the rapidity windows 0.4 < y < 0.6. We note that, in this rapidity domain, the thermal 
contribution is practically zero. So, we can conclude that the proposed model satisfactorily describes the m±- 
spectra even in those regions, where the contribution from a pure thermal source is absent and the spectrum is 
determined only by the rescattering of net nucleons. 

All this leaves us with the continued challenge of applying the model to other experiments and problems. 
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